Roton energy gap and spontaneous symmetry breaking 
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We study elementary excitations in superfluid helium-4 employing an approach based on the spon- 
taneous symmetry breaking. In particular, we calculate the roton energy gap at zero temperature 
A(0). The relation that we have derived is A(0)/fcs?A = 4. The theoretical value of A(0)/fes is 
8.707 K, which is significantly close to the experimental value of 8.712 K. The deviation between 
the theoretical and experimental values is less than 0.1%. 
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I. INTRODUCTION 



Recently, the author has derived the following formula 
for the A-temperature in superfluid helium-4 [lj 



(1) 



where po is the density of liquid hclium-4 at zero tempera- 
ture, m and a are the mass and the hard sphere diameter 
of a helium atom, respectively. Thus, A-temperature is 
determined only by the three physical parameters, pa,m 
and a. Among these parameters, the hard sphere di- 
ameter a is especially important, because it is a crucial 
parameter for hard-core repulsive interactions at short 
distances. Using eq. (HJ, one can obtain the theoretical 
value of T\, because the experimental values of po, m 
and a are known. Eq. ([TJ) gives T A = 2.194 K [l|, which 
is very close to the experimental value of 2.1768 K Q. 
The deviation between the theoretical and experimental 
values is less than 1%. In this meaning, our approach 
succeeded. However, in the previous paper [l[, we did 
not study the elementary excitations. 

Although the roton, which is an elementary excitation 
with finite energy gap A, has been extensively studied 
in the past, its physical nature still remains unclear [3j|. 
In this paper, we calculate the roton energy gap at zero 
temperature A(0) to understand the physical nature of 
roton. First, we calculate the roton energy gap using the 
effective potential of the mean field effective field theory. 
Second, we introduce the curvature matrix which allows 
us to calculate the energy gap directly. The relation that 
we have derived is very simple, A(0)//cb1a = 4. This 
relation gives A(0)/fcs = 8.707 K, which is very close to 
the experimental value of 8.712 K. Finally, we comment 
on an open problem about the sound speed of phonon. In 
the next section, we begin with a brief review of ref. [l[ 
for self-consistency of this paper. 



II. SPONTANEOUS SYMMETRY BREAKING 
AND A-TEMPERATURE 

The superfluid phase of liquid helium-4 is described by 
the Gross-Pitaevskii (GP) effective field theory, in which 
the cooperative state of many interacting helium atoms 
is represented with a single complex scalar field ip. GP 
theory is defined by 



(2) 



where tp is a one-component complex scalar field with 
both amplitude and phase, m is the mass of a helium 
atom, p is the chemical potential, A = 2irh 2 a/m is the 
coupling constant for hard-core repulsive interactions at 
short distances and a is the hard sphere diameter of a 
helium atom, respectively. Thus, GP theory is a non- 
relativistic version of the Goldstone model in particle 
physics. GP Lagrangian @ is invariant under the global 
U(l) transformations; ip — > e l6 ip, ip* — > e~ l6 p* with 
80/ dt = X70 = 0. Because U(l) symmetry is not gauged, 
no vector gauge field appears. 

In the GP theory, the A-transition is described as the 
spontaneous symmetry breaking of the U(l) symmetry. 
One can understand this statement as follows. The effec- 
tive potential V at the tree-level approximation is given 
by the third and the fourth terms in eq. |2j), which does 
not contain the derivatives; 



V 



X 



-pip tp- 



(3) 



As shown in ref. [l| , the A-temperature T\ is determined 
only by this potential. Here it should be noted that 
the chemical potential p becomes positive below the A- 
temperature [lj, and then the ground states become de- 
generate for T <T\. Fig. Q] shows the schematic temper- 
ature dependence of the effective potential V. 

Let's consider the ground state of this system below 
the A-temperature. In general, the ground state is given 
by solving the condition 
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dV 

dtp 



-pip* + \(<p*ip)<p = 0. 



(4) 



FIG. 1: The schematic temperature dependence of the effec- 
tive potential V. ip = is not stable for T <T\. 



FIG. 2: The potential per number density at zero tempera- 
ture. The depth represents the A-temperature. 



This condition yields two solutions, 



<p = 0, \tp\ 



(5) 



As shown in Fig.[TJ the trivial solution ip = is not stable 
for T < T\, but another solution \<p\ = y/ p/X is stable. 
Furthermore, because \(p\ 2 is equal to the number density 
n of superfluid, the following relation is satisfied below 
the A-temperature, 



(6) 



m = \lj = 



where p s is the superfluid density and m is the mass of 
a helium atom. From eq. ©, the chemical potential at 
zero temperature p(0) is given by 



m 



(7) 



where po is the total density of liquid hclium-4 at zero 
temperature, and we have used the relation p s (0) = po- 

Before deriving the formula for T\, we should see the 
spontaneous symmetry breaking. For T > T\, the effec- 
tive potential has only one minimum at <p — 0, and then 
the global U{V) symmetry is not broken. For T < T\, 
the field ip rolls down from "hill-top" at <p — to the 
bottom, and then the global U(l) symmetry is sponta- 
neously broken. This phenomenon is the A-transition. 

The derivation of the formula for T\ is very easy as fol- 
lows. When T = 0, the scalar field ip is located at the bot- 
tom, and then it takes the nonzero value \<p\ = \J p(0)/A. 
Hence the global U(l) symmetry is spontaneously broken 
in this phase. As the temperature increases above the A- 
temperature, the scalar field tp is located at tp = 0, and 
the global U(l) symmetry is recovered. Therefore, it is 
enough that we consider how energy is necessary to take 
the scalar field cp from the bottom to "hill-top". From 
eq. ([31), the depth of the potential Vq is given by 



V = -V(\p\ = v/MOVA) 



M(0)no 



(8) 



where no = Po /m is the number density of liquid hclium- 
4 at zero temperature. The depth of the effective po- 
tential represents the required energy density to recover 



the global U(l) symmetry. Because the scalar field ip has 
two degrees of freedom (its amplitude and phase), the 
A-temperature is given by 



Vo/no - -k B T x x 2. 



(9) 



Fig. [2] shows this situation. From eqs. ©, ®, © and 
A = 2irh 2 o~ /m, we obtain the following formula for the 
A-temperature, 



T x = 



p Trh 2 a 



(10) 



Using eq. (fit)]) , one can calculate the A-temperature, 
because the experimental values of po? m an d o~ are 
known. Substituting the following values into eq. (E 



Po 
rn 
a 



0.1451 g/cm 3 , 
6.6465 x 10" 24 
2.639 A, 



one obtains the theoretical value 



T\ = 2.194 K. 



(11) 
(12) 

(13) 



(14) 



This is very close to the experimental values of 2.1768 K, 
and the deviation between the theoretical and experimen- 
tal values is less than 1 %, Thus, spontaneous symmetry 
breaking is a key concept to understand the A-transition. 



III. ROTON ENERGY GAP 

In this section, we calculate the roton energy gap 
at zero temperature A(0) using the effective potential ([3]). 

(i) linear decomposition 

Because we are interested in the excitations near the 
ground state, it is convenient to define 



tp(x,t) = ip c + T)(x,t) +ix(x,t), ip c = Jj, (15) 

where n and x are rea ^ scalar fields. This decomposition 
is known as linear decomposition. Plugging (|15p into l[5]). 



we obtain the effective potential in terms of rj and \ 

V=-£+ 2 W 2 + 2^V(^ 2 + X 2 )V +^(V 2 + X 2 ) 2 -(16) 

In eq. (|16|) . the first term (at zero temperature) deter- 
mines the A-temperature as shown in section 2, 



Vo = 2 ^ = k B T x n . 



(17) 



Using the relation /«(0)/A = n (see, eq. dU)), we obtain 

M (0) = 2k B T x . (18) 

This is the direct relation between the chemical potential 
at zero temperature and the A-temperature. 

However, for the present purpose, the second term in 
eq. (flU)) is more important, because the coefficient of if 
term represents the potential energy of the real scalar 
field rj. Therefore, the coefficient of rj 2 should be identi- 
fied with the energy gap of roton at zero temperature 

A(0) = 2/x(Q). (19) 

From eqs. (fT5)) and p^|) . the roton gap A(0) is given by 

A(0) = 2/i(0) = Ak B T x , (20) 



A(0)/k B T x = 4. 



(21) 



Substituting T\ — 2.1768 K into eq. [[21"]). one obtains 



A(0)/k B = 8.707 K 
A(0)/* B = 8.712 K 



(theory), 
(experiment). 



(22) 
(23) 



Thus, our relation (|2ip is very precise. The deviation 
between the theoretical and experimental values is less 
than 0.1%. Fig. [3] shows the roton energy gap as a func- 
tion of temperature. From fig. [3l one can find that the 
relation A{Q)/k B T\ = 4 is in good agreement with the 
experimental values at low temperatures (T < 0.8 K). 

Here we comment on the energy gap of phonon. 
Because the coefficient of x 2 term in eq. (| 16[) is zero, 
the energy gap of the scalar field \ is exactly zero. This 
means that phonon is gapless, and it is a consequence of 
the Nambu-Goldstone theorem. 

(ii) polar decomposition 

Although the linear decomposition is enough to calcu- 
late the roton energy gap, we consider another decompo- 
sition for completeness 



. (24) 



where a and tt are real scalar fields. This is known as 
polar decomposition. Because the field 7r is divided by ip c , 
the fields a and 7r have the same physical dimensions. In 




1 1.5 
temperature [K] 

FIG. 3: The roton energy gap as a function of temperature. 
The theoretical value is in good agreement with data. 



this decomposition, the field ir is the Nambu-Goldstone 
mode and a is the Higgs mode. 

Plugging ([24]) into ([3]) , we obtain the effective potential 
in terms of a 



2v /vt 3 + \j\ 



(25) 



Thus, we find that the Nambu-Goldstone mode 7r van- 
ishes in the effective potential. The field tt appears only 
in the kinetic terms in the Lagrangian. Therefore, the 
Nambu-Goldstone mode or phonon is clearly gapless. 

The physical meaning of the 1st and 2nd terms in 
eq. (f2"B"|) is the same as that of linear decomposition. The 
1st term determines T\, and the 2nd term determines 
A(0). Thus, the linear decomposition and the polar de- 
composition yield the same physical results; phonon is 
gapless and the roton gap is given by A(0)/k B T\ = 4. 



IV. CURVATURE MATRIX 

In this section, we introduce the curvature matrix, 
which allows us to calculate the energy gap of elemen- 
tary excitations more directly. 

First, let's consider the effective potential at T = 



u(o) = - A1 (o)^> + ^(^v) 2 . 

Here we introduce the 2x2 curvature matrix C 

d 2 V(0) 



(26) 



C 



I d 2 V(0) 

V 



d<p*d<p 
d 2 V(0) 



d<p*dtp 



<fi=v>Ao) 



¥>=Vc(0) 



d<p*dtp 
d 2 V(0) 



<P=<Pc(fl)j 



(27) 



where <p c (0) is the field value at T = (see, eq. (J5J), 



M (0) 



A 



(28) 



Because the matrix C represents the curvature at the 
bottom of the effective potential, two eigenvalues of the 



4 




wave number k (A 1 ) 



FIG. 4: Phonon-roton dispersion curve. 

matrix C should be identified with the energy gap of 
elementary excitations at zero temperature. 

From eq. (|26p . the curvature matrix C is given by 

C = MO) (\ = 2k B T x (] . (29) 

Therefore, two eigenvalues of the matrix C are 

0, Ak B T x . (30) 

These eigenvalues represent the energy gap of elementary 
excitations; 

A phonon = (Nambu-Goldstone mode) (31) 
A roton = 4k B T x (Higgsmode). (32) 

Thus, using the curvature matrix, one can calculate the 
energy gap of elementary excitations very easily. 



where 777,4 is the mass of a helium-4 atom and /i(0) is the 
chemical potential at T = 0. Although the original Bo- 
goliubov theory cannot calculate the sound speed quanti- 
tatively, we can calculate it because we have derived the 
relation (fl~8|) . From eqs. (fl~8]) and (|33]) . one obtains 

y 7774 V "7,4 

where we have used 7774 = 6.64645 x 10~ 24 g. This value 
is clearly inconsistent with the experimental value of 238 
m/s. Thus, although the original Bogoliubov theory ex- 
plains the linear dispersion of phonon (e = cp = chk), it 
fails to explain the sound speed quantitatively. 

Here we propose a possible expression of sound speed, 
which is consistent with the experimental value. In 
eq. (|33|) (of Bogoliubov theory) , if we replace 7774 by the 
effective roton mass m* = 7774/271 = 0.16 7774, one obtains 

c=M=,/^= 238 m/, (35) 
y 777, V m* 

Surprisingly, this is completely consistent with the exper- 
imental value as shown in fig. 01 but its physical origin 
still remains to be solved. This indicates that the ro- 
ton effective mass (or equivalently, the roton kinetic en- 
ergy) should be considered in order to calculate the sound 
speed of phonon. However, at present, the theoretical ex- 
planation for eq. ([33)1 is not known. We hope that this 
open problem will be discussed in the community. 



PHONON 



VI. CONCLUSION 



Finally, we comment on an open problem about the 
sound speed of phonon. 

According to the standard Bogoliubov theory Q , it is 
well known that the sound speed of phonon c is given by 



'Mo) 
7774 



(33) 



In conclusion, we have calculated the roton energy gap 
at zero temperature using an approach based on the 
spontaneous symmetry breaking. The relation that we 
have derived, A(0) jk B T x — 4, is in good agreement with 
the experimental value. We hope that this work con- 
tributes to the progress of low temperature physics. 
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